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ABSTRACT 

We report the detection of a deep broad H I 21 cm absorption system at z = 0.78 toward the radio 
source [HB89] 2351+456 (4C+45.51) at z = 1.992. The H I absorption was identified in a blind spectral 
line survey conducted at the Green Bank Telescope spanning 0.63 < z < 1.10 toward a large sample of 
radio continuum sources. The H I column density is Nhi — 2.35 x 10 19 (T S //) cm -2 , where T s is the spin 
temperature and / is the continuum covering factor of the absorbing gas. For T s / f > 8.5 K, this system 
is by definition a damped Lya absorption system (Njji > 2 x 10 20 cm -2 ). The line is unusually broad, 
with a FWHM of 53 km s _1 and a full span of 163 km s _1 , suggesting a physically extended H I gas 
structure. Radio surveys identify damped Lya systems in a manner that bypasses many of the selection 
effects present in optical/UV surveys, including dust extinction and the atmospheric cutoff for z < 1.65. 
The smooth broad profile of this H I 21 cm absorption system is similar to the z — 0.89 H I absorption 
toward PKS 1830—211, which suggests that the absorber toward [HB89] 2351+456 is also a gravitational 
lens and a molecular absorption system. However, very long baseline interferometry and Hubble Space 
Telescope observations show little evidence for gravitational lensing, and BIMA millimeter observations 
show no HC0 + (1 — > 2) or HCN (1 — > 2) absorption down to r = 0.15 (3er) in 5 km s _1 channels. 
Although this radio damped Lya selection technique would include dusty, molecule-rich systems, [HB89] 
2351+456 appears to be a "vanilla" H I 21 cm absorber. 

Subject headings: galaxies:ISM — quasars:absorption lines — quasars: individual [HB89] 2351+456) 



1. INTRODUCTION 

Identification of damped Lya (DLA) systems, quasar 
absorption line systems with neutral hydrogen column 
density Nhi > 2 x 10 20 cm~ 2 , has historically been an 
optical/UV pursuit. DLA catalogs suffer from the selec- 
tion effects associated with an optically selected sample 
of redshifted Lya lines, including the atmospheric cutoff 
at z < 1.65, the requirement of a strong UV continuum, 
and the potential loss of dusty DLA systems. To overcome 
the UV atmospheric cutoff, low ionization absorption lines 
such as Mg II A2796, 2803 and Fe II A2600 provide prox- 
ies for DLA systems that can subsequently be observed in 
the Lya line from space or in the H I 21 cm line from the 
ground (Rao & Turnshek 2000; Lane 2000). The H I 21 cm 
line, unlike Lya, does not saturate, is unaffected by dust, 
and can occasionally offer multiple lines of sight through 
the absorption system. H I 21 cm absorption is a prerequi- 
site for molecular absorption, so a dust-independent DLA 
selection method maximizes the likelihood of adding to the 
very small sample of extragalactic molecular absorption 
systems. These systems provide insight into the physical 
state of the cold interstellar medium in galaxies and can 
constrain the cosmic evolution of physical constants (e.g., 



Wiklind & Combes 1997; Darling 2004; Kanekar, Chen- 
galur, & Ghosh 2004). While there does not appear to be 
a dominant undetected population of high column density 
DLA systems compared to the DLA population overall, 
optically selected samples may exclude a large fraction of 
the dusty absorption line systems that do exist (Carilli et 
al. 1998; Ellison et al. 2001), thereby significantly dimin- 
ishing the likelihood of identifying molecular absorption 
systems. 

Traditional approaches to H I 21 cm observations of 
DLA systems have focused on identifying the radio line 
toward optically identified systems of known redshift (e.g., 
Carilli et al. 1996), although there have also been "blind" 
narrow instantaneous bandwidth surveys for intervening 
H I absorption and surveys for intrinsic absorption in 
quasars (e.g., Brown & Roberts 1973; Carilli et al. 1998). 
New radio astronomy facilities equipped with broadband 
receivers and high resolution spectrometers now offer the 
opportunity to conduct line searches for DLA systems in 
a manner similar to traditional optical/UV DLA surveys, 
spanning large redshift intervals with high spectral resolu- 
tion. For example, the Green Bank Telescope 1 (GBT) of- 
fers continuous frequency coverage from 290 to 1420 MHz 
[z = 0-3.9 in H I), large instantaneous redshift coverage, 
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and high velocity resolution (e.g., Az = 0.5 and Av = 5 
km s- 1 at 800 MHz). 

We have conducted a DLA survey toward strong ra- 
dio sources at the GBT at 800 MHz spanning the redshift 
range z = 0.63-1.10 with 5 km s -1 resolution that is sensi- 
tive to all DLA systems. We select sources with predicted 
780 MHz flux density > 1.7 Jy as interpolated from 1.4 
GHz White & Becker (1992) and 365 MHz Douglas et 
al. (1996) measurements. This is not a perfect radio sur- 
vey: we also require sources to have z > 1.06, which biases 
the sample toward systems with optical identifications and 
measurable lines. 

We report the first detection of the Green Bank DLA 
survey, a deep broad H I absorption system at z = 
0.78 toward the radio source [HB89] 2351+456 (hereafter 
2351+456, but also known as 4C+45.51, B3 2351+456, 
NVSS J235421+455304, and WMAP 074, among oth- 
ers). To investigate the nature of this absorption system, 
we present BIMA millimeter observations of HCO + and 
HCN. We also reexamine published optical spectra and 
infrared photometry and find that the H I absorption sys- 
tem provides new interpretations of observed properties of 
2351+456. This Letter assumes the concordance cosmol- 
ogy: H = 71 km s" 1 Mpc" 1 , Cl M = 0.27, and fi A = 0.73. 

2. OBSERVATIONS 
2.1. Green Bank Telescope 

We observed 2351+456 at 773.75 MHz spanning 192.5 
MHz at the GBT on 2004 January 26. We observed four 50 
MHz IFs with 2.5 MHz overlap in two linear polarizations, 
spanning the range 677.5-870 MHz or z HI = 0.63-1.10. 
Each bandpass was broken into 8192 6.1 kHz channels, 
but adjacent channels in the autocorrelation spectrometer 
are not independent and the effective spectral resolution 
is 12.2 kHz or 4.7 km s _1 at the band center (z = 0.84). 
Fast sampling (1.5 s data records) was used to facilitate 
radio frequency interference (RFI) excision. Observations 
of 2351+456 include a single 5 minute position-switched 
observation and a 25 minute on-source scan. 

An H I absorption line was identified at 798.2 MHz 
(z = 0.78) in individual 1.5 s records during observations. 
We find no evidence for RFI contamination of the detected 
21 cm line, and we find no evidence that the line is an 
RFI signal in the off-source scan or an out-of-band har- 
monic; the HI absorption line is clearly visible in individ- 
ual records, and this line is not seen in any other program 
source. 

Records were individually calibrated and bandpasses 
flattened using the winking calibration diode. The 
25 minute on-source scan was calibrated in a position- 
switched mode using the 5 minute off-source scan at the 
cost of added noise. A 25 minute scan on [HB89] 2338+042 
observed during the same session as 2351+456 was cali- 
brated in an identical manner as a template for subtraction 
of the feed resonance from the 2351+456 spectrum (Fig- 
ure 1, inset). A least-squares fit of the template spectrum 
to the 2351+456 spectrum in line- and RFI- free regions 
indicated that a small shift (0.073 Jy) and scaling (0.975) 
was adequate to fit the resonance in 2351+456. We sub- 
tracted the shifted and scaled template from the 2351+456 

2 The AIPS++ (Astronomical Information Processing System) is fr 
may be obtained from http://aips2.nrao.edu 



spectrum to obtain a resonance and baseline subtracted 
spectrum. A polynomial baseline was subtracted from a 
fit of 5 MHz centered on the H I line to remove bandpass 
ripples of order 2 MHz (the total span plotted in Figure 
1 is 1.6 MHz). We reached a typical rms noise of 4 mJy 
in 6.1 kHz (2.3 km s^ 1 ) channels. All data reduction was 
performed in AIPS++. 2 

2.2. BIMA HCO+ and HCN Observations 

We observed 2351+456 in the HCO+ J = 1 -> 2 
(178.3751 GHz) and HCN J = 1 -> 2 (177.2611 GHz) lines 
redshifted to 100.2344 and 99.60678 GHz, respectively (as- 
suming z = 0.7796), at the Berkeley-Illinois-Maryland As- 
sociation (BIMA; Welch et al. 1996) array at Hat Creek 
on 2004 April 9 and 18. We used correlator mode 6, which 
provides four independent spectral windows. Two of them 
are 100 MHz wide with a spectral resolution ~ 1.56 MHz 
(~ 5 km s^ 1 ), and two are 200 MHz wide with a resolu- 
tion ~ 3.12 MHz (~ 10 km s _1 ). The two high-resolution 
windows were partially overlapped and centered on each 
expected absorption feature, covering a velocity range of 
Av <~ ±200 km s _1 around it. The two low-resolution 
windows extended this velocity coverage to Av ~ ±750 
km s _1 . 

Observations of a passband calibrator were interleaved 
(a few minutes every hour) with the science observations 
and used to remove passband features. During the first 
track on April 9, we used QSO 0102+584, which was con- 
siderably weaker (S u ~ 1.5 Jy) than predicted by the flux 
database and introduced some noise into the passband cor- 
rection. For the second track on April 18, we switched to 
QSO 0136+478, which yielded better results (S u ~ 2.8 
Jy). The passband-corrected spectra, which are feature- 
less, are shown in Figure 2. The sensitivities achieved after 
passband correction are a ~ 35 mJy in the high resolution 
region. Like most flat spectrum sources, 2351+456 ap- 
pears to be variable, with a flux density of S v w 719 mJy 
on April 9 and S v « 648 mJy on April 18. 

3. RESULTS AND ANALYSIS 

Confusion, coupled with RFI-induced effects like slopes 
and steps in the spectral bandpass, indicates that the GBT 
spectral baseline does not accurately reflect the contin- 
uum emission from 2351+456. Within the 15.5' (FWHM) 
GBT primary beam at 800 MHz, there is a contribution of 
roughly 100 mJy at 1.4 GHz from sources of unknown spec- 
tral index (Condon et al. 1998). We thus interpolate the 
800 MHz flux density from previous interferometric mea- 
surements: 1.873+0.056 Jy at 1.4GHz (Condon et al. 1998) 
and 2.158 ± 0.041 Jy at 365 MHz (Douglas et al. 1996) 
indicate a flat-spectrum source with a = —0.11 ± 0.03 
[S v oc v a ). Hence, at 798.2257 MHz, the expected (and 
henceforth assumed) flux density is 1.987 + 0.066 Jy in the 
absence of variability (but see, e.g., Aller, Aller, & Hughes 
(2003)). 

Fits to the H I profile (Figure 1) require a minimum of 
seven Gaussian components to obtain residuals consistent 
with the spectral noise. These fits are not presented be- 
cause they are not unique and thus not physical. The seven 
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component requirement and the width of the absorption 
system indicates that the H I profile is a highly blended, 
perhaps continuous, ensemble of gas complexes sampling 
a galaxy-scale range of velocities and physical dimensions. 

The peak optical depth in the H I line is 0.320 (545.0 
mJy) at z = 0.779452 ±0.000014. We use the most conser- 
vative error estimate of the redshift of ±1 channel in lieu of 
unique reliable Gaussian fits to the profile. The integrated 
optical depth is 13.05 km s , and the column density is 
thus N HI = 1.8 x 10 18 (T S //) frdv = 2.35 x 10 19 (T S //) 
cm" 2 , where T s is the spin temperature in K and / is the 
covering factor of the continuum source by the absorbing 
gas. The absorption system qualifies as a DLA system 
(Nhi > 2 x 10 20 cm" 2 ) for T s /f > 8.5 K. Typical val- 
ues span T s = 100-1000 K (Kanekar & Chengalur 2003), 
indicating a column density from Nhi — 2.35 x 10 21 to 
2.35 x 10 22 cm" 2 for / = 1. The FWHM of the absorption 
line is 53 km s _1 , and the total span of absorption is 163 
km s _1 . Since the line is unusually broad and the width 
is due to bulk velocities rather than to thermal kinetic 
broadening, the usual constraint on the kinetic tempera- 
ture T fe < (1.2119 x 10 2 FWHM 2 )/(81n2) w 61,000 K is 
not physical. 

In the BIMA spectra, no HCO+ (1 -> 2) or HCN 
(1 — > 2) lines were detected with |r| > 0.15 (105 mJy) 
at 3(J significance in 5 km s" 1 channels (Figure 2). In con- 
trast, the four known molecular absorption systems show 
large optical depths, from r = 0.32 in HCO + (1 — > 2) 
and t = 0.25 in HCN (1 -> 2) toward B3 1504+377 to 
saturated lines toward PKS 1830-211 and B0218+357 
(Wiklind & Combes 1995, 1996a,b, 1997). The limit 
on HCO + absorption toward 2351+456 corresponds to 
N H co+ < 3.2 x 10 13 cm" 2 , which translates into the limit 
Nh 2 < 1-6 x 10 22 cm" 2 assuming a gas kinetic tempera- 
ture of 50 K and N H2 = 5 x 10 8 N H co+ (Liszt & Lucas 
2000). For T fe = 25 K, we obtain N H2 < 5 x 10 21 cm" 2 , 
which is equivalent to Ay = 2.5 for a Galactic dust-to-gas 
ratio. 

4. DISCUSSION 

The nature of the absorption system toward 2351+456 
remains unclear. The H I 21 cm absorption profile shows 
striking similarities in width and lack of prominent sub- 
structure to the absorption toward PKS 1830—211 (Chen- 
galur, de Bruyn, & Narasimha 1999). The absorber toward 
PKS 1830—211 is both a gravitational lens and a molecular 
absorber (Subrahmanyan et al. 1990; Wiklind & Combes 
1996a), and the unusually broad and smooth H I profile 
likely represents illumination of a significant fraction of an 
H I spiral disk by a background radio source magnified by 
gravitational lensing (Patnaik, Muxlow, & Jauncey 1993; 
Chengalur, de Bruyn, & Narasimha 1999). But BIMA ob- 
servations show that the 2351+456 absorber does not have 
strong molecular absorption lines, and existing VLBI and 
Hubble Space Telescope (HST) observations can explore 
the possibility of gravitational lensing of 2351+456. 

VLBI observations of 2351+456 at 18 cm (1664 MHz) 
show a ~ 20-30 mas resolved core plus a one-sided jet 
morphology with total extent ~ 110 mas (Polatidis et al. 
1995). There is no clear evidence of a second image on 
scales of tens to ~ 150 mas in the VLBI map (although 
there is a tantalizingly strong knot of emission ~ 110 mas 



from the core that could be a weak counterimage) , and an 
archival /-band (F814W) HST WFPC2 image shows no 
second image on scales > 200 mas. Although scales rel- 
evant to gravitational lensing are probed by HST images 
and VLBI maps from ~ 10 mas to arcseconds, no clear 
companion images have been identified. We conclude that 
2351+456 is unlikely to be gravitationally lensed. 

The physical scale subtended by the 18 cm VLBI core 
at the absorber redshift is 0.15-0.22 kpc, but this scale is 
likely to be larger at 37.5 cm (800 MHz). The absorption 
line width strongly suggests that the H I absorption arises 
from a unusually large physical scale, of order the size 
of the illuminating "beam" of the VLBI core. The high 
frequency picture of 2351+456 is dramatically different: 
Lister (2001) shows that 2351+456 at 43 GHz subtends 
only ~ 1 mas. The millimeter illuminating beam is thus 
only ~ 8 pc at the absorber redshift, which can potentially 
explain the lack of molecular absorption in the BIMA spec- 
tra. Although the scales probed at centimeter wavelengths 
are large, the millimeter scales are in the "pencil beam" 
regime that makes the detection of molecular absorption 
so rare and difficult. 

Previous optical observations of 2351+456 may offer 
some clues to the nature of the Hi absorption system at 
z = 0.78. Based on the 1987 September 15 optical spec- 
trum of 2351+456, which shows a dramatic continuum in- 
crease from 1986 October 26 but preserves flux in emission 
lines and the equivalent width in the Mg II absorption line 
at z = 1.99 (Lawrence et al. 1996), we identify additional 
absorption lines consistent with two absorbing complexes. 
One is roughly at the redshift of the quasar (z w 1.99), 
showing possible absorption in Lya, Fe II AA2344, 2374, 
2382, 2586, and 2600, Mg II AA2796 and 2803, and Mg I 
A2852 (many of these are blended). The other is consis- 
tent with the redshift of the H I detection at z = 0.78, with 
absorption in Fe II AA2344 and 2600, Mg II AA2796 and 
2803, and possibly H7 and H<5 (many of these are blended, 
and H<5 roughly coincides with Fe II A2586 at z — 1.99, 
so one is likely a spurious identification). Some of these 
absorption lines can also be seen in the Stickcl & Kiihr 
(1993) spectrum. Rao & Turnshek (2000) show that half 
of absorption systems with Mg II A2796 and Fe II A2600 
equivalent widths greater than 0.5 A are DLA systems and 
that low redshift DLA systems that do not meet these re- 
quirements are rare (only one is known). Hence, although 
a higher resolution spectrum is required to extract reliable 
absorption redshifts and equivalent widths, the presence 
of Mg II and Fe II absorption lines is consistent with the 
identification of a DLA system at z = 0.78. 

The absorption system toward 2351+456 may also have 
been mistakenly identified as the host galaxy of the quasar 
itself. Kukula et al. (2001) modeled the HST KLCMOS H- 
band surface brightness of 2351+456 and obtained a half 
light radius of 18.5 kpc (corrected to the concordance cos- 
mology) for the quasar host assuming an r 1 / 4 model (de 
Vaucouleurs 1948). If the extended emission is in fact from 
the intervening absorber at z — 0.78, then the half light 
radius of the absorber is 16.2 kpc. If this is the case, then 
a de Vaucouleurs profile may not be an appropriate pa- 
rameterized fit. 

5. CONCLUSIONS 
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We have identified a DLA system toward 2351+456 via 
the H I 21 cm absorption line'. We identify Mg II AA2796 
and 2803 and Fe II AA2344 and 2600 lines in the optical 
spectrum of Lawrence et al. (1996), in support of the H I 
21 cm absorption line at z = 0.78 being a DLA system. 
Although the H I absorption line profile of 2351+456 is re- 
markably similar to the profile of PKS 1830-211, the DLA 
system toward 2351+456 does not appear to be either a 
molecular absorption system or a gravitational lens. This 
new detection provides a proof-of-concept for blind H I 
surveys with broadband radio receivers and backends that 
offer a new dust- and atmosphere-independent method to 
identify DLA systems. Broadband radio surveys for DLA 
systems are perhaps key to identifying exceptionally rare 
molecular absorption systems and may identify gravita- 



tional lenses as well. 
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Fig. 1. — The H I absorption spectrum toward 2351+456. Zero velocity refers to barycentric redshift z = 0.779452, and the velocity scale is 
in the rest frame of the absorption system. Zero flux density and zero optical depth correspond to the continuum level of 2351+456, assumed 
to be 1.99 Jy. Inset: Position-switched spectra of 2351+456 (bottom) and [HB89] 2338+042 (top) centered on the 800 MHz receiver feed 
resonance at 796.8 MHz. The spectrum of [HB89] 2338+042 has been shifted by +0.5 Jy. [HB89] 2338+042 provides a template for the feed 
resonance and general bandpass subtraction of the 2351+456 spectrum. The narrow features on the red side of the resonance are RFI. 
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Fig. 2. — Millimeter-wave spectra of 2351+456 obtained at BIMA for the redshifted transitions of HCO+ (1 — > 2) and HCN (1 -> 2) toward 
2351+456. The thin trace corresponds to the spectral region observed at high resolution, while the thick gray trace corresponds to the low 
resolution correlator windows. The thick black line at the bottom illustrates the expected Vlsr of the molecular absorption for z = 0.779452, 
taking the heliocentric to LSR frame corrections into account. 



